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Abstract 

We study the decay width and forward-backward asymmetry of the lepton pair for 
the inclusive decay b — > sr+r^ in the two Higgs doublet model with three level flavor 
changing neutral currents (model III) and analyse the dependencies of these quantities 
on the model III parameters, including the leading order QCD corrections. We found 
that there is a considerable enhancement in the decay width and neutral Higgs effects are 
detectable for large values of the parameter . 
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1 Introduction 



Currently, there is an impressive experimental effort for studying rare B-meson decays at SLAC 
(BaBar), KEK (BELLE), B- Factories, DESY (HERA-B) since these decays are rich phenomeno- 
logically. They are induced by flavor changing neutral currents (FCNC) at loop level in the 
Standard model (SM) and with the forthcoming experiments, it would be possible to test the 
flavour sector of the SM in a high precision, as well as to reveal the physics beyond, such as 
two Higgs Doublet model (2HDM), Minimal Supersymmetric extension of the SM (MSSM) [||], 
etc. 

Among the rare B decays, B K*l^l^ process has received a great interest since the SM 
prediction for its branching ratio (Br) is large enough to be measured in the near future. This 
decay is induced hj b sl^l~ transition at the quark level and in the literature it has been 
investigated extensively for / = e,/i in the SM, 2HDM and MSSM g]- |T5|. When / = e,/i, 
the neutral Higgs boson (NHB) effects are safely neglected in the 2HDM because they enter in 
the expressions with the factor me(^)/mvi/- However, for I = r, this factor is not negligible and 



NHB effects can give important contribution. In |T^, |T^, B — > X^r+r^ process was studied in 
the 2HDM and it was shown that NHB effects are sizable for large values of tanf3. 

In this work, we study the h —>■ st^t~ decay in the general 2HDM , so-called model HI. 
We include NHB effects and make the full calculation using the on-shell renormalization pre- 
scription. We investigate the dependencies of the differential decay width dT/ds and the decay 
width r on the scale invariant lepton mass square " s" and some model HI parameters, namely 
mH±, ^^jjj, and Further, we calculate the differential (direct) forward-backward asymme- 

try Afb{s) {Afb) of the lepton pair in terms of the above parameters. We show that a large 
enhancement is possible in the decay width of the process b —>■ st~^t~ for some values of the 
model HI parameters and NHB effects become considerable for large values of ^^^t-. 

The paper is organized as follows: In Section 2, we present the leading order (LO) QCD 
corrected effective Hamiltonian and the corresponding matrix element for the inclusive b —>■ 
ST^T^ decay. Further, we give the expression for Apsis) and Aps of the lepton pair. Section 
3 is devoted to the analysis of the new Wilson coefficients coming from the NHB effects and 
the dependencies of dT/ds, F, Afb{s) and ApB on the the Yukawa couplings ^^^65 ^N,TTy the 
charged Higgs mass , the parameter s and to the discussion of our results. In Appendices, 
we give the explicit forms of the operators appearing in the effective Hamiltonian and the 
corresponding Wilson coefficients. 
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2 The inclusive b sr^r decay in the model III 

Model III (2HDM) permits the flavour changing neutral currents in the tree level and the prize 
is various new parameters, i.e. Yukawa couplings. These couplings are responsible for the 
interaction of quarks and leptons with gauge bosons, namely, the Yukawa interaction and in 
this general case it reads as 



Cy = V^jQiL<l)lUjR + r]j'QiL<f)iD,R + &QiL(t)2U,R + ^^jQiL<i)2DjR + h.c 



(1) 



where L and R denote chiral projections L{R) = 1/2(1 =f 75), 0^, for k = 1,2, are the two 
scalar doublets, QiL are quark and lepton doublets, Ujji, Djji are the corresponding singlets, 
Tlfj^ , and ^^j^ are the matrices of the Yukawa couplings. The Flavor changing (FC) part of 
the interaction is given by 



^Y,FC = ifiQiL4>2UjR + ^^QiL<p2DjR + h.c. . 



(2) 



The choice of 0i and 02 



V2 







+ 



( V2x' 



with the vacuum expectation values. 



J_( V2H+ 
v/2 Iv Hi + tH2 



(3) 



(4) 



ensures decoupling of the SM and beyond. In eq.(|^) the couplings ^^'^ for the FC charged 
interactions are 



^ch — ^neutral VcKM , 



^ch 



VcKM ^neutral , 



where Cneutrai Q is defined by the expression 



cU,D 



U,DYl^U,DyU,D 



(5) 



(6) 



Here the charged couplings appear as linear combinations of neutral couplings multiplied by 
VcKM matrix elements (see |]T8|] for details). 

Now we would like to start with the calculation of the matrix element for the inclusive 
h ST^T^ decay. The procedure is to integrate out the heavy degrees of freedom, namely 



^In all next discussion we denote ^, 



neutral 



as ^ 



U,D 
N ■ 
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t quark, ,H^,H^,Hi, and H2 bosons in the present case and obtain the effective theory. 
Here denote charged, if". Hi and H2 denote neutral Higgs bosons. Note that Hi and H2 
are the same as the mass eigenstates and in the model III respectively, due to the choice 
given by eq. (^). The QCD corrections are done through matching the full theory with the 
effective low energy one at the high scale /x = mw and evaluating the Wilson coefficients from 
my/ down to the lower scale /i ~ 0{mf,). In the model III (similar to the models I and II, 
2HDM) neutral Higgs particles bring new contributions to the matrix element of the process 
b ST^T^ (see eq. (|23|) ) since they enter in the expressions with the mass of r lepton or related 
Yukawa coupling As being different from the model I and II, in the model HI, there exist 

additional operators which are the flipped chirality partners of the former ones. However, the 
effects of the latter are negligible since the corresponding Wilson coefficients are small due to 
the discussion given in section 3. Therefore, the effective Hamiltonian relevant for the process 
b —>■ ST^T~ is 



n 



eff 



'^W''^'' |EC.(/i)0.(/i) +Y.CQMQ^M , (7) 



where Oj are current-current {i = 1,2), penguin (z = 3, ...,6), magnetic penguin (i = 7,8) and 
semileptonic {i = 9, 10) operators. Here, Ci{fi) are Wilson coefficients normalized at the scale 
fi and given in Appendix B. The additional operators Qi{i = 1,..,10) are due to the NHB 
exchange diagrams and Cq^{^) are their Wilson coefficients (see Appendices A and B) . 

During the calculations of NHB contributions, we use the on-shell renormalization scheme 
to overcome the logarithmic divergences. Taking the vertex function 

yRen / 2n ^ -pO / 2\ , yC fo\ 
neutr \f I ^ neutr \f J i ^ neutr i V / 

and using the renormalization condition 

^neutriP ~ ^neutr) ~ 0' (9) 

we get the counter terms and then calculate r^|"^^(p^) . Here the phrase neutr denotes the 
neutral Higgs bosons, H^, and and p is the momentum transfer. 

Now we give the QCD corrected amplitude for the inclusive b —>■ st^t~ decay in the model 

HI, 

M = ^V^ifel^,:|c9^^^(s7;.PL5)r7'^r + Cio(s7^PL&)r7'^75r 

- 2C'7^(s«a^.j9'^PK6)r7^r + C'Q,(sP«6)fr + C'Q,(sP^6)f75r| . (10) 



Using Eq.(|T^), the differential decay rate reads as 

ds ^ U7c^f{mc/mby ' \ s ) ^ ' ^ ' 

witti 



l + 2s + — (l-4s) 
s 



+ l2Re{CrCl"*) il + —j + ^\CQf{s - At') + ^\Cq,\' + 6Re{C^oC*Q,), (12) 

wfiere s = p'^/ml, t = rrir/mb, and /(x) is tlie pliase-space factor given by /(x) = 1 — + 
— x® — 24x^ logx. In tlie above expression for the differential decay rate, we use the inclusive 
one since, in the heavy quark effective theory, the leading terms of inclusive decay rates of the 
heavy hadrons in the l/mt expansion becomes that of the free heavy quark, 6-quark in our 
context. 

The forward-backward asymmetry ApB of the lepton pair is another physical quantity which 
can be observed in the experiments and provide important clues to test the theoretical models 
used. Using the definition of differential Aps 

A r \ Jo dsdz J-i dsdz /i o\ 

JO "-^dsdz + J-1 "-^dsdz 

with z = cos 6, where 6 is the angle between the momentum of the b-quark and that of r"*" in 
the center of mass frame of the dileptons r+r^, we get 

AM^) = (14) 

Here, 

E{s) = Re{Cl^^Cl,s + 2CrCl, + Cl"C*Q^t + 2CrC*Q^t) (15) 
In addition, ApB can be defined as 

ApB = ^ —. (16) 

Note that during the calculations of F and ApB, we take into account only the second 
resonance for the LD effects coming from the reaction b sipi st^t^, where i = 1,..,6 
and divide the integration region for s into two parts : < s < ('"'^2 g^j^j {m^,^+o.o2) ^ 

b b b 

s < 1, where m^^ = 3.686 Gel^ is the mass of the second resonance (see Appendix B for LD 
contributions) . 



3 Discussion 



In the general 2HDM model, there are many free parameters, such as masses of charged and 
neutral Higgs bosons and the complex Yukawa couplings, i^j^ ■, where i, j are quark flavor indices 
and these parameters should be restricted using the experimental measurements. Usually, the 
stronger restrictions to the new couplings are obtained from the analysis of the AF = 2 (here 
F = K, Bdi D) decays, the p parameter and the B decay. 

The neutral Higgs bosons and Aq give contributions to the Wilson coefficient C7 (see the 



appendix of [T^ for details) 



C^"{mw) = {Vt,V*)-' E ^NM^N,s^^, (17) 

where rrii and Qi are the masses and charges of the down quarks (z = d, s, h) respectively. These 
expressions show that the neutral Higgs bosons can give a large contribution to the coefficient 
Cj which is in contradiction with the CLEO data [ppj], 

Br{B X,7) = (3.15 ± 0.35 ± 0.32) 10"^ . (18) 

Such dangerous terms can be removed by assuming that the couplings ^^^^ {i = d,s,b) and 
^N,db small enough to be able to reach the conditions ^N,bb^N,is ^ ^N,db^N,ds ^■ 
The discussion given above results in the following restrictions: ~ and ^^^j ~ 0, where 



the indices i,j denote d and s quarks . Further using the constraints [21], coming from the 
AF = 2 mixing, the p parameter and the measurement by CLEO Collaboration eq. (0) 



we get the condition for ^Ntc, ^Ntc << ^nu ^^^^ ^^^'^ account only the Yukawa couplings 



^N,tt 

of quarks ^^j^ and ^^bb- As for we do not consider any constraint and increase this 

parameter to enhance the effects of neutral Higgs boson. (For further discussion about the 



restrictions of the model HI parameters see [18, 21 



In this section, we study the Wilson coefficients CQ^{mb) and CQ^{mb) coming from NHB 
effects and s, and ^^.^t- dependencies of dT/ds and F for the inclusive decay b st^t^ , 
restricting [C7'^''^[ in the region 0.257 < [C7'^'^[ < 0.439 due to the CLEO measurement, eq.(|T^) 
21|[ for details). Our numerical calculations based on this restriction and throughout these 



see 



calculations, we use the redefinition 



we take the scale p = rrib and use the input values given in Table (HD. 
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Parameter 


Value 


rrir 


l.Yo (GeVj 


rric 


1.4 (tjeVj 


nib 


4.0 (^(oreV j 


"em 


1 on 


At 


U.U4 




U.iUo it U.Ui 


rat 


i ( (^Vore V j 




ou.zu ^^vjrt; V J 


mz 


91.19 (GeV) 


^QCD 


0.225 (GeV) 


as{mz) 


0.117 


sinOw 


0.2325 



Table 1: The values of the input parameters used in the numerical calculations. 

In Fig. 15 (0), we present rriho dependence of CQ,{mb) for C^^^ > 0, I'^bh = 40 m,, (Sm;,), 
^Ntt = ^GeV in the case \rtb\ = < 1 {nb > !)• Here Cn^imb) lies in the region bounded 

by solid lines. For rtb > 1, fnho = 80 GeV and m//o = 100 GeV^, the value of CQ^{mb) changes 
between 0.020 and 0.045. However for rtb > 1, we get values, -9 and -12, more than two orders 
of magnitude larger compared to ones for \rtb\ < 1, for the same value of rriho. Since CQ^(mfe) is 
directly proportional to ^^^t, its value may further increase with the increasing values of ^^tt- 



The corresponding 2HDM model II value of Cq-^ (nib) can be extracted from as beeing ~ 0.4 
for large tan /3, tan /3 = 25. 

For completeness, in Figs. |]and we give rriHO dependence of Cq-^ {rrib) and m^o dependence 
of Cq^imb), for Cj^^ > 0, ^§bb = 40 m;,, ^^,tt = 5GeV in the case \rtb\ < 1. As seen from 
Fig. ^, niAo dependence of CQ^{mb) is relatively weaker and for niAo = 80 GeV, Gq^imb) is 
between nearly 0.0284 and 0.0291. For rtb > 1, G^^^ > 0, ^§^bb = Sm^ and ^^ ..^ = 5 GeV, 
Gq^imb) reaches up to the value of —0.38. The 2HDM model II value of \Gq^{mb)\ is ~ 0.4 for 
tan/3 = 25 [0]. 

Now we continue the analysis of the measurable quantities F and Aps oi the process under 
consideration. In the following, we use the numerical values m/^o = 150 GeV , niho = 80 GeV 
and m^o = 80 GeV in our calculations. 

In Fig. 1^, we plot the differential F of the decay b — > sr+r^ with respect to the parameter s 
for = 40 mb, = I GeV and charged Higgs mass mH± = 400 GeV in case of the ratio 
\rtb\ < 1. Here the differential F lies in the region bounded by dashed (small dashed) curves for 
G^^-^ > {Gy^^ < 0). A small enhancement is possible especially for G^^^ > case compared to 
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the SM (solid curve). Further, the restriction region of the differential F for model III becomes 
narrower with increasing or decreasing values of the parameter s. Fig. ^ is devoted the same 
dependence of the differential F including the long distance (LD) effects. Here Cj^'^ < case 
for model III almost coincides with the SM (solid curve). In case of the ratio > 1, extremely 
large enhancement, 3 orders larger compared the \rtb\ < 1 case, is reached even for the small 
values of ^^ ^^ (see Fig. 0). 

I'D 

Fig. I shows dependence of F of the decay under consideration for C,^ .^^ = 1 GeV and 
charged Higgs mass mH± = 400 GeV in case of the ratio \rtfy\ < 1. Here F is almost non-sensitive 
to However for rtt > 1 case (Fig. ^ F is strongly sensitive to for Cj > 0. Further, 
F is 2 orders (3 orders) larger compared to the SM result for C^-^'^ < {Cy"^-^ > 0) even for 

''N.bb ^ 2 

Fig. |TD| is devoted to the dependence of F to the charged Higgs mass mj^. F has a weak 
dependence (almost no dependence) on for C^^^ >0 {Cj" <0). 



For completeness, in Figs. |TT] and |T2| we also present C^rr dependence of F for large 
values of ^^t,-. Sensitivity of F to ^^.n- increases with the increasing values of this parameter. 
F enhances for extremely large values of ^^tt ^^"^ ^^is is the contribution due to the NHB 
effects. For \rtb\ < 1 the NHB effects are small and destructive up to the large values of ^^^t, 
^^^^^ = 800 GeV. For Gj^^ > 0, ^ = 40 and ^^^^^ = 1(100) GeV this effect is at the order 
of the magnitude %0.1 (4) of the overall contribution. However, it is positive for rtb > 1 and it 
becomes considerable with increasing values of ^§tt- -^^^ C7 > 0, the small value = 3 
and ^^..^ = 1(100, 200) GeV, the NHB contribution can reach the magnitude %0.15 (7, 26) of 
the overall contribution. 

Our results on Apb{s) and Aps for the decay under consideration are presented through 
the graphs given by Figs. |T3]- |T3|. In Fig. ^ Apb{s) is shown for ^§f,b = 40 m^, ^^^^ = 1 GeV 
and charged Higgs mass mH± = 400 Ge^ in case of the ratio \rtb\ < 1. Here Apb{s) lies in 
the region bounded by solid lines for Gj^^ > 0. Dashed line presents Gj^^ < case and the 
SM result coincides with this line. There is possible negative values of Apb{s) due to the LD 
effects. For rtb > 1, Apb{s) almost vanishes (~ 10~^). 

Fig. Ill is devoted to dependence of ApB for C,MrT = I GeV, charged Higgs mass 

I— I ' 



= 4:00 GeV and \rtb\ < 1. Here, ApB is not sensitive to especially for large values 



of this parameter. The SM and model HI average resuhs for Gj" < {G'/^ > 0) are 0.340 
and 0.340 (0.325), respectively. ApB is sensitive to the parameter for its small values in 
the case where rtb > 1 and G^^^ < (Fig. |l^). The enhancement over the SM is possible for 
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< 0.4, namely Aps can reach the value of 0.45. The restriction region for Afb is large for 
this case. However, for Cj-^-^ > 0, Afb almost vanishes. 

The NHB effects on Afb is sensitive to the couphng ^^.^^ as it should be. For \rtb\ < 1 and 
Cy^^ > 0, the NHB contribution is -% 0.15 for ^^^^ = 1 GeV and % 1.2 for ^^.^^ = 100 GeV in 
case the parameter = 40. Increasing ^^,tt causes to the enhancement in the NHB effects. 
For m > 1 and Cj^^ < 0, the NHB effects are negative and it increases the overall result by 
% 10 for ^ = 0.4 and ^^^^^ = 100 GeV. For rtb > 1 and C^^^ > 0, the NHB effects to Afb 
are negligible. 

Now, we would like to summarize our results. 

• r for the process under consideration is at the order of 10~^ for \rtb\ < 1 and C^-^^ > 
results is greater compared to Cj-^-^ < one. On the otherhand, for > 1, there is 
a considerable enhancement, three order larger compared to the SM case even for small 
values of Further, F is not sensitive to for \rtb\ < 1, however strong sensitivity 
to this parameter is observed for rtb > 1- 

• ApB is not so much sensitive to the model 111 parameters for \rtb\ < 1- For r^f, > 1 , 
there is a possible enhancement in the Afb for small values of however it becomes 
negligible with increasing 

• The NHB effects becomes important for the large values of the Yukawa couphng C^ ,-^. 

Therefore, the experimental investigation of F and Afb ensure a crucial test for new physics 
and also the sign of C^^^ . 
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Appendix 



A The operator basis 

The operator basis in the 2HDM (model III ) for our process is [l^, |2^, ^ 



02 = {sLa1^lCLa){cLpl^hL|3), 

03 = {sLali^bLa) J2 i^LfsYgLp), 

q=u,d,s,c,b 

q=u,d,s,c,b 
q=u,d,s,c,b 

Oq = {sLa1t,hLp) iqRfdl^QRa), 
q=u,d,s,c,b 

g 

g 

g 

q=u,d,s,c,b 
q=u,d,s,c,b 
q=u,d,s,c,b 

Qe = Yf^(^l&^) E i^RQl) 

q=u,d,s,c,b 

Q7 = ^2(''1^"'^r) E 

q=u,d,s,c,b 

Qs = ^2(''l^"^bR) E (4^, A) 

q=u,d,s,c,b 



Q9 = ^2(''L^"'b^R) E ifL^.uq'^R) 

q=u,d,s,c,b 

QlO = ^2(''L^"'b^R) E i4^,^QL) (19) 

q=u,d,s,c,b 

where a and (3 are SU (3) colour indices and JF^*^ and Q'^'^ are the field strength tensors of the 
electromagnetic and strong interactions, respectively. Note that there are also fiipped chirality 
partners of these operators, which can be obtained by interchanging L and R in the basis 
given above in model III. However, we do not present them here since corresponding Wilson 
coefficients are negligible. 

B The Initial values of the Wilson coefficients. 



The initial values of the Wilson coefficients for the relevant process in the SM are 

C'JU^w) = 0, 

C^'im^) = 1, 
ySM / \ _ — 2Xj — 8xf — 5xf + 7Xt 



A{xt - 1)4 ^ ' 2A{xt - 1; 



3 



5 ■> 



^« ^""^^ - ~4(x,-l)4^^'^*+ 8(0;, -1)3 ' 

sin'^Uw sm 9 

C^ir(^H^) = -4^iB{xt)-C{x,)) , 
sm b'vK 

Cgfl^w') = 2 = 1,. .,10. (20) 
The initial values for the additional part due to charged Higgs bosons are 

Ctsimw) = 0, 

C^imw) = Y'F,{yt) + XY F^iyt) , 

C^{mw) = Y^Gi{yt) + XY G2{yt) , 

C^imw) = Y'HM , 

C^oi^w) = Y^L^{yt) , (21) 



where 



~ ^ (^N,bb ^N,sbTr ) 

nih \ ' VtbJ 



rrib \ ' Vtby 
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and due to the neutral Higgs bosons are 



32'K^m%mtQi{zA, ytj^bivt) 



MO 



'S>i{zA,yt)Q5{yt) 



). 



^Q2 (^N,tt) 



-xyxtyt 



{x - l)xt{yt/zA - 1) - (1 + x)yt 



xyxt + 2za _ 2in[ ^^®5(a:i) 



Ql{zA,Xt) 

x{yt + xt{l - yt/zA)) - '2yt 



(06 -{x- y){xt - yt)){'d3{zA) + {x- y){xt - yt)zA) 



r^A^fcD \ _ 3 ^N,rT^N,bb x^yt + 2|/(a: - l)xtyt - za{,xI + Be) , xljl - yt/zA) , ^ , f zaQ 



QsizA) 



'Q2{zAy^' 



9\^N,ttf^bmr f xtil - 2y)yt , {-1 + 2 cos^ 9w)i-l + x + y)yt 



+ 



zh{Oi{zh: yt)xyt + cos^ 9w {-2x'^{-l + xt)yy1 + xxtyy^ - Q8^g)) \ ^23^ 



^Qi {^N,tt) — 



-l + 2cos''ew)yt/ 1 



647r^m^()^t \ cos^ 9w 
'S>5iyt)zH 



o^iyt) 



67 



xtyt ^ xtyt{xy - Zh) 



Obiyt) Qi{zH,yt) 



-2xt In 



Ol{zH:yt) 



g mbiTirXt 



- 1 + 



-l + 2x)xt ^2xt{-l + i2 + xt)y) 



e^ixt) + y{l - xt) 



Acos^ ew{-l + X + y) + xt{x + y) Xt{x{xt{y - 2zh) - 4:Zh) + 2zh) 



+ 



cos^ OwO^ixt) Qi{zH,xt) 
yt((-l + x)xtZH + cos^ 9w{{3x - y)zH + Xt{2y{x - 1) - zh{2 - 3x - y)))) 



cos2 9w(Q3(zh) + x(xt - yt)zH) 



+2 {xt In 



Qb{xt)zH 
Ql{zH,Xt) 



+ In 



xjyt - xt)zH - Qz{zh) 
{Qb{xt) - xt)ytZH 



) . 
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^N,rr^N,bbi^N,ttf ( {Q^{yt)zh{yt " i)(x + y - 1) - e,{zh, yt){Q,{yt) + yt) 



2 In 



327r2m^o 
ZhQ^jyt) 



Qi{zH)Q,{yt) 



9^^N,rr^N,ttmhXt ( 2(-l + (2 + Xt)y) Xt{x - l){yt - Zh) 



+ 



x{xt{y - 2zh) - 4:Zh) + 2zh 



Q'2iZh) 

(1 + x)ytZh 



+ 2 In 



ZhQ^jXt) 
QliZh,Xt) 



Qi{zh,xt) xyxtyt + Zh{{x - y){xt - yt) - ©e) 

_^ ©9 + ytZhijx - y){xt - yt) - ©6)(2x - 1) ^ x{ytZh + Xt{zh - yt)) - 'iytZh ^^ 



^Qi {^N,bb) 



Zh&eCde - (x - y){xt - yt)) 



9 ^N,TT^N,bb ( yxtyt{xxf{yt - Zh) + QQZh{x - 2)) 



Q2{zh) 



ZhQ2{Zh)Q% 



+ 2 In 



ZhQt 



Q2{Zh) 



where 



ei(cj,A 

Q2{U0 
©2(^ 
©3(^ 

e5(A 
ee 
67 
eg 



= —{—l+y — y\)uj — x{y\ + uj — uj\) 

= {xt + y{l - xt))ytuj - xxtiyyt + {jjt - 

= Q2{'^,xt yt) 

= {xt{-l +y)- y)yti^ + xxtiyyt + + yt)) 

= 1 — X + xu 

= X + A(l — x) 

= {xt + yil-xt))yt + xxtil-yt) 

= iyiyt-l)-yt)zH + x{yyt + iyt-l)zH) 

= yt{2x\l + xt){yt - 1) + xt{y{l - yt) + yt) + x{2{l -y + yt) 

+ xt{l-2y{l-yt)-3yt))) 

= -Xt{-1 + X + y){-yt + x{2yt - \)){yt - Zh) - XtytZh{x{l + 2x) - 2y) 

+ y'ii.xt{x^ - y{l - x)) + {l + x){x- y)zh) 



(24) 



and 



Xt 



mr 



yt 



m: 



Zh 



mr 



Zh 



mr 



ml 



za 



mr 



m 



The exphcit forms of the functions F^2){yt), Gi(2){yt), Hi{yt) and Li{yt) in eg. (PT|) are given as 

yt{7~5yt~8y^) , yK^yt - 2) 



Fiiyt) 



72{yt - 1)3 + 12{y, - 1)^ 
12 



F2{yt) 
Hiivt] 



yti^yyt - 3) yt{-'iyt + 2) 

12(yi-l)2+ Q{yt-lf 

yt{-y^t + %t + 2) 



+ 



-1/t 



24(i/< - 1)3 
- 3 ) , i/t 

4(yt - 1: 



4(2/* - 1) 



+ 



2(2/t - 1)= 



1 — Asin^Ow xyt 



yt 



472/2 - 792/i 



yt 

38 



1 
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^i(2/*) 



108(2/i - 1)= 



xyt 



sin?6w 



3yf ~6yt + 4: 
myt - 1)4 
1 



+ 



(yt - ly 



(25) 



Finally, the initial values of the coefficients in the model III are 



dx 



l-x 



+ C'qi ((^iV,tt)'') + i^N,tt) + i^N,bb)) 



c, 



2HDM , 



dx 



l-x 



dy{cim 



2HDM , 



2HDM , 



mw 



mw 





rrib 






[mw) 


1 r^2HDM 
+ *-^Q2 


m-r 


sin^ 






rrib 






[mw) 


r<2HDM 


m-r 


sin^ 


Ow 


= 0, 


i = 


5,. 


.,10. 







(26) 



Here, we present Cg^ and Cq^ in terms of the Feynmann parameters x and y since the inte- 
grated results are extremely large. Using these initial values, we can calculate the coefficients 
Cf^^^\fi) and CQ^'^*^(/i) at any lower scale in the effective theory with five quarks, namely 
u, c, d, s, b similar to the SM case |jT3|, |16|, |19|, ^ |. 



The Wilson coefficients playing the essential role in this process are 



C^Q^^^fi), Cq^^^lfi) and Cq^^^^fi). For completeness, in the following we give their explicit 
expressions. 



Cr^if-^) = Cr^'^ili) + {Cr^'\l^) + iVc C6^™'(/i)) , 
where the LO QCD corrected Wilson coefficient Cj^"'^^^'^^ (fi) is given by 

^LO,2HDM^^^ = r^^^/^'C'/'^^'imw) + (8/3)(r/i4/23 - r]''^/'')C','''''\mw) 



(27) 



mw) ^iV"' 



{2i 



i=l 



13 



and r] = as{mw)/ a sifi), hi and Oj are the numbers which appear during the evaluation . 

Cg-^-^ (fi) contains a perturbative part and a part coming from LD effects due to conversion 
of the real cc into lepton pair r+r": 



(29) 



where 



cr(/z) 



and 



Y 



ri2HDM 



+ h{z, s) (3Ci(/i) + C2{^^) + 3C3(/i) + C^{^j) + SCsl/i) + ^(/i)) 

- s) (4C3(^) + 4C4(/x) + 3C5(^) + C6(/i)) (30) 

- ^/i(0, s) (C3(/x) + 3Q(/i)) + ^ (3C3(/x) + C,{^i) + 3C5(/x) + CM) , 



r ' r 



em 

(3Ci(/x) + C2(/x) + 3C3(/i) + C4(/i) + 3C,{fi) + CM) 



(31) 



In eq.(^), the functions h{u, s) are given by 

h{u, s) = 



m, 8 ^84 

m m M H — X 

27 9 



9 9 
2 



9 



(2 + x)|l -X 



1/2 



(in 



— ivr ) , for a; = 



2 arctan ■ 



< 1 



-1' 



for X = ^ > 1, 



MO,.) = ^-9l---9l-^+9^- 



(32) 



(33) 



with M = The phenomenological parameter k in eq. (|3T| ) is taken as 2.3. In eqs. (30) and 
(PI|), the contributions of the coefficients Ci(/i), Cq{ii) are due to the operator mixing. 
Finally, the Wilson coefficients Cq^{h) and Cq^{h) are given by 



(34) 
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Figure 2: Same as Fig.|T], but for ^^f,^ = 3 m;, and rtb > 1. 
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Figure 3: CQ-^{mb) as a function of m^/o for — 40 m^, Cn,tt — ^GeV, mn^ — 400 Gel^ and 
m/jo = SOGeV^ in case of the ratio \rtb\ < 1- 
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Figure 4: Cq^irnb) as a function of m^o for = 40 m^, ^^.^^ = 5GeV and mH± — 400 Ge\^ 
in case of the ratio \rtb\ < 1- 
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Figure 5: Differential F as a function of s for = 40 nib, ^n,tt ~ ^ GeV and mH± = 400 GeV 
in case of the ratio \rtb\ < 1. 
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Figure 7: The same as Fig ||, but for 



3 and rth > 1- 
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Figure 8: F as a function of 
\rtb\ < 1. 
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Figure 9: The same as Fig. R| but for rj^ > 1. 
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Figure 10: F function of mH± for ^^^^ 



40 mb, 



1 GeV in case of the ratio \rtb\ < 1- 
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Figure 11: F as a function of C,^ 
\rtb\ < 1- 
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Figure 12: The same as Fig. |Tl| for = 0.1 mb in case of the ratio > 1. 
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Figure 13: Differential ApB as a function of s for ^^^j, = AOuif,, ^^^^ — IGeV and mH± 
AOOGeV including LD effects in case of the ratio \rtb\ < 1- 
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Figure 14: ApB as a function of ^ for ^^^^ = 1 GeV, mH± = 400 GeV and \rth\ < 1- 
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